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We propose two hypotheses which characterize the quasi-equilibrium state that realizes after a 
cold collapse of self-gravitating system. The first hypothesis is the linear temperature-mass (TM) 
relation, which yields a characteristic non-Gaussian velocity distribution. The second is the local 
virial (LV) relation, which, combining the linear TM relation, determines a unique mass density 
profile as p(r) = por~ i e~ roi ' r . Although this density profile is unphysical in the central region, the 
region is just inner a few percent around the center of a bound region in cumulative mass, which 
is beyond the resolution of our numerical simulations. Hence posing two hypotheses is compatible 
to the numerical simulations for almost the whole region of the virialized bound state. Actually, 
except for this inner part, this density profile fits well to the data of cold collapse simulations. 
Two families of spherical and isotropic models, polytropes and King models, are examined from a 
view point of these two hypotheses. We found that the LV relation imposes a strong constraint on 
these models: only polytropes with index n ~ 5 such as Plummer's model are compatible with the 
numerical results characterized by the two hypotheses among these models. King models with the 
concentration parameter c ~ 2 violate the LV relation while they are consistent with the R 1 ^ 4 law 
for the surface brightness. Hence the above characteristics can serve as a guideline to build up the 
models for the bound state after a cold collapse, besides the conventional criteria concerning the 
asymptotic behavior. 

PACS numbers: 



I. INTRODUCTION 

Collisionless self-gravitating systems (SGS) eventually 
settle down to a quasi-equilibrium state through the 
phase mixing and the violent relaxation processes under 
the potential oscillation lj. This quasi-equilibrium state 
is a prototype of the astronomical objects such as galax- 
ies and clusters of galaxies. These objects often show 
universal profiles in various aspects. 

It would be true that the energy distribution function 
and the mass density distribution function are quite ap- 
propriate to characterize SGS, and actually many studies 
have been done from this viewpoint. For example, the 
density profile is found to be p oc r -4 in spherical and 
cold collapses, i.e., collapses with small initial virial ra- 
tio 0, 01 and p oc r~ 3 — r~ 4 in the cluster-pair merging 
processes 0. 

On the other hand, the velocity distribution function 
has not fully been examined so far despite its impor- 
tance , anisotropy of velocity dispersion or lower mo- 
ments of velocity distributions with respect to the de- 
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1 So far, full velocity distributions have been mainly analyzed only 
for the line-of-sight velocity profiles because they can be directly 
compared with the observational data 



viations from Gaussian distributions have been mainly 
examined (e.g. 0,0)- However, the full velocity distri- 
bution function plays an important role in identifying the 
quasi-equilibrium state, which should be properly char- 
acterized by the whole phase-space probability functions. 
Merrall et al. numerically showed that the velocity distri- 
butions in the central regions of the bound particles be- 
come Gaussian, when the bound particles are well relaxed 
0. However, since SGS are manifestly non- additive sys- 
tems the ordinary Gaussian distribution would not 
be expected for the whole bound particles. Moreover, 
such non-Gaussian velocity distributions would be firmly 
connected with the density profile characteristic to SGS 

Recently, we have simulated self-gravitating N-body 
systems by utilizing the leap-frog symplectic integrator 
on GRAPE-5, a special-purpose computer designed to 
accelerate N-body simulations [lOj. The simulations were 
performed mainly in the case with cold collapse initial 
conditions. It is certain that the processes may be too 
simplified to deduce the ultimate nature of actual SGSs, 
in which merging or tidal effects take place. However, one 
of our main goals is to establish the universal properties 
of SGS through the violent relaxation process, i.e. the 
relaxation through the oscillation of gravitational poten- 
tial. This process is rather common in actual SGSs such 
as galaxies and dark matters. In this context, the cold 
collapse is considered to be a typical case which causes 
the violent relaxation. Hence it is meaningful to examine 
the cold collapse as a first step. 
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When the system experiences violent gravitational pro- 
cesses such as a cold collapse and a cluster-pair colli- 
sion, we obtained a universal velocity distribution pro- 
file expressed as the democratic (=equally weighted) su- 
perposition of Gaussian distributions of various temper- 
atures(DT distribution hereafter) in which the local tem- 
perature T(r) is defined using the local velocity variance 
(v 2 ) as T(r) := m(v 2 )(r)/3kB, where m is the parti- 
cle mass If]. Moreover, we have found that the lo- 
cally defined temperature linearly falls down in the in- 
termediate cluster region outside the central part, pro- 
vided it is described against the cumulative mass M r , 
i.e., dT/dM r = const. This fact is consistent with the 
appearance of DT distributions. 

Furthermore to the linear TM relation, we have also 
obtained another peculiar fact that the LV relation be- 
tween the locally defined potential energy and kinetic 
energy holds except for the weak fluctuations. More 
precisely, the local temperature T(r) is proportional to 
the local potential $(r), with constant proportionality 
6fcflT(r) = — m$(r) for a wide class of cold collapse 
simulations. It was originally pointed out by Edding- 
ton that Plummer's model satisfies this condition exactly 
[T^ |. Moreover it was also pointed out that Plummer's 
model can be extended to the anisotropic model of a stel- 
lar system satisfing the LV relation [l3|. We found that 
the bound region, i.e., the region composed of the bound 
particles keeps this relation quite well at least during the 
time interval much longer than the initial crossing time. 
This is a prominent correlation in the phase space dis- 
tribution function, and also constructs another backbone 
character of SGS after a cold collapse. 

In sec Inland IIIII we will show that the above two hy- 
potheses for SGS are supported by the results of a va- 
riety of cold collapse simulations. Then in sec II VI com- 
bining the two hypotheses, i.e., the linear TM relation 
and the LV relation, we obtain a unique density profile 
p(r) — po( r Ao) which asymptotically behaves 

as p oc r~ 4 . Justification of this density profile will also 
be discussed. In sec[V] we will ascertain that the above 
two hypotheses can arise even in the collisionless quasi- 
equilibrium state, where the region inside a half-mass 
radius is not yet under the effect of two-body relaxation. 
In sec I VII the above two hypotheses will be explored in 
two typical collisionless static models with spherical con- 
figuration and isotropic velocity dispersions; polytropcs 
and King models. Sec lVIII is devoted to the summaries 
and conclusions of this paper. 



II. LINEAR TM RELATION 

Let us now further consider the linear TM relation. 
In our recent work we have examined the velocity 
distribution function of collisionless SGS after a cold col- 
lapse in numerical N-body calculations. According to this 
work, immediately after a cold collapse, well-relaxed and 
almost spherical-symmetric bound state is formed. In or- 



der to examine the local property, we use the cumulative 
mass M r within the radius r, i.e., M r := Att J q dr'r' 2 p(r') 
as a measure of the radial distance from the cluster cen- 
ter, which is defined as the position of the potential min- 
imum. In our numerical simulation, we divide the whole 
system into several concentric shells with equal width in 
the coordinate M r and consider the averaged quantities 
within each shell as local variables. 

In our simulation, particles in the outermost region, 
which includes 20% of the total mass, escape almost 
freely just after the cold collapse and are irrelevant to 
the quasi-equilibrium state. Therefore we concentrate on 
the local temperature within the inner region initially in- 
cluding 80% of the total mass. More precisely, we define 
the bound particles as the particles inside the M r , at 
which the local temperature T(M r ) takes the minimum 
value. Hereafter, we pay our attention to the bound re- 
gion composed of these bound particles. In the central 
part of this region, the two-body relaxation proceeds and 
the particles evaporate toward outside. Moreover, in this 
central region, the effect of cutoff e proposed in numeri- 
cal simulations cannot be ignored completely. Discarding 
this part from the bound region, we found, that the lo- 
cal velocity variance, i.e., the local temperature linearly 
decreases to zero in the coordinate M r , 



(v 2 )(r) = —2- = c(M tot (t) - M r ), (1) 



where c is a positive constant and M tot (t) is the total 
mass of the bound particles at t. M to t{t) is initially equal 
to 80% of the total mass but is gradually diminished, as 
several particles evaporate from the bound state. 

Here we examined the deviation of TM relation from 
the linear relation (JJJ for a variety of simulations of cold 
collapse. For the estimate of the deviation, we first take 
the time-average for the velocity dispersion at each shell 
to derive the time-averaged (u 2 )-mass relation. Then we 
linearly fit it as Q to derive (v 2 ) fa and subtract it from 
(v 2 ). In Fig^ the deviation normalized by (v 2 ) at the 
half mass radius are depicted for several simulations with 
total particle's number TV = 5000. The deviations are 
very small and are less than 10 percent at each M r for 
all of the simulations we examined. Hence, this result 
is quite robust against initial conditions with N = 5000 
provided the collapse is cold. Even for the case of cluster- 
pair collision, this linear relation is widely obtained. We 
also examined the N dependence of the linear TM re- 
lation for the case with initially homogeneous spherical 
sphere (FigJ2J|. As N increases, the particles collapse 
more and more coherently, which causes a shock region 
moving slowly outward after a collapse. Although such a 
shock region enhances the deviation from the linear TM 
relation, the region within a shock region satisfies the 
linearity quite well. 
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FIG. 1: A mass dependence of velocity dispersion 
({v 2 (M r /Mtot))) obtained by a typical cold collapse simula- 
tion; spherical collapse (SC) and cluster-pair collisions (CC). 
In the case of SC, 5000 particles are distributed with a power 
law density profile (p oc r~ a ) within a sphere of radius 7? and 
the initial virial ratio (bt) is set to be small. In the case of 
CC, each cluster has the equal number of particles (2500) and 
all particles are homogeneously distributed within a sphere of 
radius R and is set to be virialized initially. The initial sepa- 
ration of the pair is 6R along the x axis. In all of the simula- 
tions, softening length e = 2~ S R is introduced to reduce the 
numerical error caused by close encounters. The inset figure 
shows that locally averaged velocity dispersion is plotted as 
a function of M r (t) / M to t{t) for SC(&t = 0, a = 0) as a typical 
example. Apparent linearity is remarkable. The broken line 
is the best fit line with c = 9.8G/R (Eq.Q). The open circles 
represent the time averaged numerical data from t = 5tff and 
t = 100t//, where tff is the initial free fall time defined as 
tff := yj R 3 /GM to t(0)- The deviation from the best fit line 
(< v 2 >fu) normalized by the velocity dispersion at the half 
mass is shown for some numerical simulations with different 
initial conditions. The value of each deviation is under 10%. 



III. 



LOCAL VIRIAL RELATION 



Besides the linear TM relation, there is another char- 
acteristic for SGS as we now argue. It is well known that 
the gravitationally bound system approaches a virialized 
state satisfying the condition 



W + 2K = 0, 



(2) 



where W and K are, respectively, the averaged potential 
energy and kinetic energy of the whole bound system. 
This is a global relation that holds for the entire system, 
after the initial coherent motion fades out. 

Here we define the locally averaged potential energy 
and kinetic energy inside the radius r, respectively as 



W r 



K r ■■ = 



1 



$(r')p(r')47rr' 2 dr', 



r V)(r') 



p(r')4:TTr' dr' , 



(3) 



(4) 



FIG. 2: A mass dependence of velocity dispersion 
({v 2 (M r j Mtot))) obtained by a typical cold collapse simula- 
tion SC (bt = 0, a = 0) with different particle number N. 
The case with N = 5 x 10 3 , 10 4 , 5 x 10 4 , 10 5 are depicted. 
The peak appears in the case with TV = 5 x 10 4 , 10 5 around 
M r ~ 0.7M to t, which represents a shock region. 



1.5 



0.5 



I 

SC( b,=0, a=0) 
SC ( b t =0.1, a=0 ) 
SC( b t =0.3,a=0) 
SC( b t =0.5, a=0) 



i-B-i 



SC( b t =0, a=l) 
SC( b t =0, a=2) 
CC 




0.2 



0.4 0.6 

M r / Mtot 



0.8 



FIG. 3: The LV relation for some numerical simulations with 
N=5000. The local virial ratio b is plotted as a function of 
M r j Mtot- The initial condition of each simulation is the same 
as FigQ The virial ratios of each shell are time averaged from 
t = 5t// and t = 100t//. 



where *(r') means the local object 7k- evaluated at r'. 
Then we extend the virial relation @ locally as 



W r 



2K r = 0, 



(5) 



and examine how precisely this relation is locally attained 
inside a bound state. 

More precisely, the above relation Eq.JSJ) is equivalent 
to the purely local relation 

2(v 2 )(r) = -$(r), (6) 

at each position r. Hence, we define the LV ratio 

b(r) := -2(« 2 )(r)/$(r), (7) 
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FIG. 4: The LV relation for some numerical simulations with 
different particle number N. The local virial ratio b is plotted 
as a function of M r /M to t- The case with N = 5 x 10 3 , 10 4 ,5 x 
10 4 ,10 5 are depicted. The initial condition of each simulation 
is the same as Fig|5] The virial ratios of each shell are time 
averaged from t = 5t// and t = 100t//. The peak appears in 
the case with N = 5 x 10 4 , 10 s around M r ~ 0.7M tot , which 
represents a shock region. 



and examine the value of b(r) for each shell. For a wide 
class of collapses including cluster-pair collision examined 
in the value b(r) takes almost unity: it deviates from 
unity less than 10 percent upward in the intermediate 
region and a little more downward in the inner and outer 
region for all of the simulations with N = 5000(Fig[3J). 
Hence a wide class of cold collapse simulations with N = 
5000 yield the relation ijfjj quite well. We also examined 
the N-dependence of the LV relation for the cases in FigQ] 
We got the result that the LV relation is affected by the 
shock as well as the linear TM relation in the case with 
N > 5 x 10 4 , where the LV ratio b takes the peak at the 
shock region, the value b becomes smaller than unity 
around the shock region in order to compensate for the 
enhance at the shock region. However, b oscillates around 
unity even in such cases, which supports the LV relation 
in the zeroth order. Hence, this LV relation should be 
an another characteristic of SGS, and this is our second 
hypothesis. 



IV. DENSITY PROFILE 



We have explored the two hypotheses of SGS so far, 
i.e., the linear TM relation Eq.QJ and the LV relation 
Eq.JfjJ for the bound cluster. Combining these two rela- 
tions, we now derive a mass density profile. We simply 
substitute Eq.© into Eq.lQJ to obtain 2c(M tot - M r ) = 
— $. Differentiating this with r and using Poisson's equa- 
tion, we obtain a differential equation for the mass den- 



sity p(r): 

d . a . G 2 
which admits the unique solution 
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p(r) = Po(y) e-*°' r , 



(8) 



(9) 



where r<j :— G /2c and po is a constant. The cumulative 
mass M r for this density profile is easily derived as 



M r = Mt nt e- ro /\ 



where M to t is the total mass and is described as 



M tot = 47rrgp . 



(10) 



(11) 



The central part of this density is unphysical, since it 
increases with increasing radius r inside r* := tq/A. This 
fact warns that the two hypotheses we proposed are not 
always justified in the entire region. However, M r inside 
e r* is just 



M r (r*) = e' 4 M tot - 0.0183M tot , 



(12) 



from IjlOl) . Hence this unphysical region is at most 2 
percent of the full bound region in cumulative mass. 

Even if we put a constant density at this inner unphys- 
ical region and continuously connect it to the outer ana- 
lytical solution JHJ, this inner region is at most three per- 
cent of the full bound region in mass (See Appendi jXj l . 
which is inner than the smallest value of M r for numer- 
ical data points in Fig^and Fig03 Hence we can safely 
neglect this unphysical region. Actually the function © 
fits well to our numerical results (Fig[5J). The asymptotic 
behavior p[r) cx r -4 is consistent with previous works 
[l4l Il5| . In sec lVIl we will see how this central part is 
modified with several collisionlcss static solutions. 



V. COLLISIONLESS OR COLLISIONAL? THE 
PROPERTIES OF THE SYSTEM COMPATIBLE 
WITH THE TWO HYPOTHESES 

In the previous section, we numerically showed that 
the bound system after a cold collapse retains the two 
hypotheses during a long time interval on average. Here 
a natural question arises, i.e., do these two character- 
istics occur just after a cold collapse or are they grad- 
ually developed through the two-body relaxation? Now 
we check this issue by estimating the two-body relaxation 
time scale t re i, which is described as 



-3/2 



0.065 LJ- 



mG*(p) log(l/e) 



(13) 



where (v 2 ) and (p) are the velocity variance and the den- 
sity averaged inside the half mass radius, respectively 
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FIG. 5: A log-log plot of a typical density profile after a cold 
collapse. The initial condition of this simulation is the same 
as SC (bt = 0, a = 0) in Fig0 The open circles represent the 
numerical data at t — 100*,,. The solid line is the density 
profile (Eq.Q) with the best fit parameter (Log 10 (ro/J?) = 
— 1.25), which corresponds to M r \ r =r / 'Mtot ~ 0.43. 



FIG. 6: A mass dependence of velocity dispersion (< 
v 2 (M r ) >) for the case of a spherical collapse (run SC(N — 
5000, bt = 0.5, a — 0)), for which ** ei is almost 41*,,. The ve- 
locity dispersion of each shell at the different time (* = 5*,,, 
10*,,, 50tff, and 100*,,) are superposed. 



TABLE I: The relaxation time for a bound state t* el is de- 
picted. The initial condition of each simulation is a spherical 
collapse (SC). Each parameter N, bt, and a denote the par- 
ticle number, the initial virial ratio, and the initial density 
profile p ~ r~ a respectively. 
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[Tfil |. Although t re i is a function of time *, the value set- 
tles down to a constant value t* el when the system stays 
in a quasi-equilibrium bound state. From the definition 
of t re i(t), the system inside a half mass radius naively 
attains two-body relaxation at t > t* el . In Table[I] we 
show the value of t* el for several simulations. 

For example, because of the high concentration of the 
bound particles just after a collapse, t* el becomes very 
short in the case with a cold collapse from a homogeneous 
sphere (run SC(iV = 5000 and 50000,67 = 0, a = 0)). 
In this case, we cannot regard the bound state after a 
collapse as a collisionless system, since t* el is so short that 
the two-body relaxation is rapidly achieved just after a 
cold collapse. On the other hand, t* el becomes longer by 
increasing the initial virial ratio or by changing the initial 
particle configuration to be inhomogcncous. Especially, 
in the case of a mild collapse with initial virial ratio 0.5 
(run SC(iV = 5000, F t = 0.5, a = 0)), t* rel becomes as long 
as 40*//. 

Here with this mild collapsing case, we examined if our 
two hypotheses are achieved even before the two-body re- 



laxation is developed in the bound state. In the case with 
the linear TM relation, the slope keeps almost constant 
even at the time in which the two-body relaxation has 
not been attained (Fig@ except for the slight shift of the 
slope in the outer part at M r ~ 0.7M tot . This constant 
slope is conserved even after the half of the bound state 
are under the effect of two-body relaxation (* > 41*//). 
Hence the character of the linear TM relation appears 
during the period that the system is almost collisionless 
and is conserved even after the two-body relaxation is 
attained in more than half of the full bound region. 

The LV relation b(M r ) = 1 is also attained even in the 
early stage, except for larger fluctuations in the outer 
part M r > 0.7 Mtot (FigEJ, although this deviation has 
nothing to do with the two-body relaxation. Since the 
half of the bound state is not relaxed at t = 10*// in 
Table|U we can claim that the LV relation is attained 
quite well even in the collisionless region. 

In summary, both the linear TM and the LV relation 
appear at the stage in which the system can be regarded 
as collisionless, although in the outer part the slope of 
TM is slightly shifted and b(M r ) fluctuates a little larger. 
The fluctuation occurs at almost the same region as the 
one where the slope slightly shifts in TM relation. Wc 
ascertained that the velocity dispersion is anisotropic at 
almost the same region in these simulations. Hence, we 
suspect that both the shift of the slope and the large 
fluctuation in the outer part are closely related to the 
anisotropy of the velocity dispersion in this region. In 
the next section, on the basis of these results we will 
examine two typical collisionless static models of SGS. 
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FIG. 7: The local virial relation for run SC(7V = 5000, b t = 
0.5,o = 0) for which t* el is almost 41t//. The ratio 
b(M r /M to t) is plotted as a function of M r . The virial ra- 
tios of each shell at the different time (t = 5tff, Wtff, 50t//, 
and WOtff) are superposed. 



VI. COMPARISON WITH SPHERICAL AND 
ISOTROPIC STATIC MODELS 

We have extracted the essence of cold collapses of SGS 
in the form of two hypotheses, i.e., the linear TM rela- 
tion Eq.Q and the LV relation Eq.©. How are they 
justified? 

A collisionless SGS can be described by the Vlasov 
equation for the phase space distribution function /(r, v) 
in the mean field limit. Especially any function of the 
integral of motion can be a stationary solution of the 
Vlasov equation. Here we treat several collisionless static 
and spherical solutions of SGS and examine which class 
of solutions meet the two hypotheses we proposed. 

In a cold collapse, the velocity dispersion of a bound 
system becomes radially anisotropic especially in the 
outer regions In the preceding section, however, 

we showed that the two hypotheses are achieved quite 
well even in the inner isotropic and collisionless region 
for a mild collapsing case. Hence, as zero-th order de- 
scription, we regard the bound state as an isotropic 
system and compare it with the spherical and static 
model with an isotropic velocity dispersion described 
as f{E), in terms of the particle energy per unit mass 
E(r, v) :— <&(r) +v 2 /2. Here we pay our attention to two 
sorts of solutions, polytropes and King models as the ex- 
amples of those models with a flat boundary condition at 
the center. We compare them with the typical numerical 
data which meet both of our two hypotheses, that is, SC 
(b t = 0,a = 0) with N = 10 4 . 

First, we examined polytropes. Polytropes are the dis- 
tribution functions with the form /(r, v) = E(r, tj)™~ 3 / 2 , 
where e is the relative energy defined as E := —E + $*, 
is chosen so as to guarantee / > for E > and n 
is the polytropic index |17| . Polytropes with n > 5 ne- 
cessitates a boundary to make the total mass finite since 



it diverges at infinity. Here, since we are interested in 
the case with no boundary, we focus our attention to the 
case with n < 5. 

For polytropes with n > 1/2, it is easy to find that the 
local temperature T(r) is connected to the local potential 
$(r) as, 

(u 2 )(r) = = 3 

m 1 + n 

which is equivalent to 

^ = - l) . (15) 

W r n + 1 \2W r ) 

(See Appendix 151 for the derivation.) Since = for 
n = 5, the LV relation Eq.J^J) or © is attained for n = 5, 
which means that the LV ratio Q becomes unity for n = 
5. This special solution with n = 5 can be analytically 
expressed as l|G4(l and is known as Plummer's model [l7j . 

In FiglHl we depicted the LV ratio b as a function of 
M r for n < 5. The data plots of numerical situation are 
located almost everywhere on the constant line of n — 5 
Plummer's model except for the outer region. As the 
value of n decreases, the slope of function b(M r ) becomes 
steeper and deviates more from the LV relation |JBJ. 

TM relation of Plummer's model is obtained as (|G8|) 
and the slope dT/dM r diverges as M r — > 0. This dis- 
crepancy is plausible, since is unphysical in the inner 
region. However the deviation from the linear TM rela- 
tion is conspicuous only in the central part of the bound 
state, i.e., M r /Mtot < 0.05. This is consistent with the 
result that the unphysical region of © is at most in- 
ner five percent of the total mass. The slope of Plum- 
mer's model becomes almost constant and the data plots 
are located almost everywhere in the middle region from 
0.2M tot < M r < 0.8M tot . In fact, the coincidence of the 
Plummer's model to to the numerical data is remark- 
able, when we normalize both of the physical variables 
with the unit of G = = M to t, where is a half-mass 
radius of the bound system. (FigEland FigEB- 

However, the polytropes with n ~ 5 fail in keeping 
the TM relation in the outer region since T(M r ) falls off 
steeply there (FigEJ). Recalling that the constant slope 
of TM relation up to the outer region is a key gradient for 
the DT distribution |ll| , the polytropes with n ~ 5 fail in 
reproducing the velocity distribution characteristic after 
a cold collapse. From (|_B3|I and l|C4|l , the mass density of 
Plummer's solution behaves as oc r -5 in the outer region. 
This is also inconsistent with the numerical results and 
with the behavior of 

Secondly, we examined King models. King models 
have often been a ppl ied to fit the photometries of ellipti- 
cal galaxies (e.g. [jJllljj)- They are conventionally ana- 
lyzed with the concentration parameter c :— log 10 (r t /r c ), 
where r t and rv are the tidal radius and the core radius, 
respectively [T3. It is shown that the King model with 
c ~ 2.25 fits the observational i? 1 / 4 law quite well [l9| . 
Here we examined the LV relation and the TM relation 
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for King models in the range 0.5 < c < 3.5. For the fam- 
ily of such a parameter range, b(M r ) becomes smaller 
in the central part for larger value of c, which means 
that the potential energy is locally more dominated than 
the kinetic energy in the central part. By contrast, for 
c ~ 0.5, the virial ratio becomes flat close to the center 
and takes the similar form to those of polytropes with 
n ~ 5fFig[TT)l. As for the TM relation for King models, 
T(M r ) of the numerical result fits well to the case with 
c ~ 0.5 in the inner part M < 0.6M tot but shifts to the 
case with c ~ 2.7 in the outer part fFig ll2[) . The den- 
sity profiles of King models with c > 2.0 also fit well to 
the numerical data in the outer part, where it behaves 
as p ~ r~ 3 — r~ 4 . This is consistent with the fact that 
the case with c ~ 2.25 fits to the observational i? 1 / 4 law 
quite well. In contrast, the density profile for the case 
with c ~ 0.5 deviates from the numerical results remark- 
ably. 

Finally, we estimate the deviations from LV relation 
and linear TM relation for polytropes and King models 
more quantitatively. Here we define the deviations from 
LV relation as S v := yf (((b — l) 2 )) where ((*)) denote 
the average of the functions of M r defined as 



M, 



tot Jo 



M tot 



dM r 



(16) 



For the definition of the deviation from linear TM rela- 
tion, we first linearly extrapolate the function T e (M r ) as 
the form of T e (M r ) = AM r + B. Then we define the 
deviation as 



T\/2 



(17) 



where T\/2 is the local temperature at the half mass ra- 
dius. In the parameter range we examined, we plotted 
both 5 V and 6 t for polytropes and King models (FigEJ. 
We also calculated the S v s derived from the numerical 
simulations of cold collapse. They are obtained by av- 
eraging (b(M r /Mtot) — I) 2 at each M r in the same way 
as the integral (|16J) . We got the result that the S v de- 
rived from the simulations in FigQ] and Fig|3] takes the 
value 0.106 ±0.053, which is depicted as the gray zone in 
FigEl The parameter range whose 5 V exceeds this gray 
zone is considered to be incompatible with the numerical 
simulations in Fig|3] 

It is obvious that the S v almost vanishes for the poly- 
tropes n ~ 5, since n = 5 Plummer's model exactly 
satisfies the LV relation, but for n becomes small, 8 V 
becomes larger and exceeds the gray zone for n < 3 in 
FigEl Hence polytropes with n < 3 cannot follow the 
LV relation as well as the results of numerical simula- 
tions. On the other hand, the StS for polytropes are al- 
most constant for any parameter range 0.5 < n < 5 and 
the value is almost 0.05. As for King models, the 5 V is 
larger than 0.1 for all of the parameter range and has a 
peak around c — 2. Especially, the S v of 1.6 > c > 2.9 
is larger than 0.2, which obviously exceeds the gray zone 
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FIG. 8: LV relation for polytropes with the unit of G = rn — 
Mtot = 1. The solid curve refers to the LV ratio b(M r /Mtot) as 
a function of cumulative mass M r /M to t for Plummer's model. 
The other curves refer to b(M r / Mtot) for polytropes with n = 
3.0 (dot-dashed), n = 1.5 (dashed), n = 0.51 (dotted). The 
plots represent the numerical data of SC (bt = 0, a = 0) with 
N = 10 4 . 



in Fig^J This is in contrast to the fact that the surface 
brightness of King models fits well to i? 1 / 4 law in this pa- 
rameter range. On the other hand, the 5tS of King models 
are less than 0.1 except for the high concentration case 
c ^ 3. The strong deviation for c <^ 3 is obvious, since 
the central temperature becomes much higher than those 
in the outer side for this parameter range. 

In summary, we can see the remarkable difference be- 
tween the polytrope sequence and the King sequence as 
the collisionless static system. In the case with King, the 
higher value of parameter 3.0 > c > 2.0 is favorable to 
the density profile or linear TM relation but is not for the 
LV relation. The lower parameter range 0.5 < c < 1.0 is 
favorable both to the linear TM relation and to the LV 
relation but is not for the density profile. Hence there is 
no parameter range in King model which is favorable to 
all of the characteristics of the bound state after a cold 
collapse. On the other hand, as for the polytrope, the 
parameter value n ~ 5 is favorable to all of those char- 
acteristics in the bound state after a collapse. Note that 
we can see the difference when we pay attention not only 
to the density profile but also to the velocity information 
such as the LV relation or linear TM relation. 

Concerning the detail of the behaviour of physical 
quantities for n ~ 5 polytropes with the centrally flat 
boundary condition such as Plummer's model, they fit 
quite well to the numerical results in the inner region, 
M r < 0.6M tot . In the outer region, however, both the 
density profile and the TM relation deviate from the nu- 
merical data. We will comment on these failures of the 
Plummer's model or n ~ 5 polytropes in the outer part 
from a viewpoint of the anisotropy of a velocity disper- 
sion in the concluding remarks. 
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FIG. 9: Temperature versus cumulative Mass (M r ) for poly- 
tropes with the unit of G — — Mtot = 1. The solid curve 
refers to the local temperature ksT/m for Plummer's model 
The other curves refer to ksT/m for polytropes with 
n — 3.0 (dot-dashed), n = 1.5 (dashed), n = 0.51 (dotted). 
The plots represent the numerical data of SC (fit = 0, a = 0) 
with AT = 10 4 . 
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FIG. 10: The density for polytropes (IC6t with the unit of 
G — ru = Mtot = 1. The solid curve refers to the density 
for Plummer's model. The other curves refer to the density 
for polytropes with n = 3.0 (dot-dashed), n = 1.5 (dashed), 
n = 0.51 (dotted). The plots represent the numerical data of 
SC (fit = 0, a = 0) with N = 10 4 . 



VII. CONCLUDING REMARKS 

In this paper, we have proposed (a) the linear TM 
relation and (b) the LV relation as the hypotheses of 
self-gravitating structures after a cold collapse. Based 
on these two hypothetical proposals, we have uniquely 
determined the mass density profile @ in the quasi- 
equilibrium state. 

Our density profile behaves as p ~ r~ 4 in the outer 
region and fits well both to the data of cold collapse sim- 
ulations. Although this density profile is unphysical in 
the inner region, where it decreases with the decreasing 
radius, we showed that such an unphysical region is just a 



b 




FIG. 11: The same as Fig[H]but for King models with the unit 
of G = r h = Mtot = 1. The solid curve refers to b(M r /M to t) 
for King model with c = 0.67. The other curves refer to 
b(M r /M t ot) for King models with c = 2.12 (dashed), c = 2.72 
(dotted). The plots represent the numerical data of SC (fit = 
0,a = 0) with N = 10 4 . 
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FIG. 12: The same as Fig|^]but for King models with the unit 
of G = rh — Mtot = 1 ■ The solid curve refers to fcsT/m for 
King model with c = 0.67. The other curves refer to fcsT/m 
for King models with c = 2.12 (dashed), c = 2.72 (dotted). 
The plots represent the numerical data of SC (fit = 0, a = 0) 
with N = 10 4 . 



few percent around the central region of the bound state. 
Hence, we can safely neglect this region at least for the 
resolution of our numerical results. 

In the central region, we also have to notice the relax- 
ation process in such a central region. In fact, we exam- 
ined how two-body relaxation is developed in a bound 
state after a cold collapse. We estimated the time scale 
of the two-body relaxation with t* el , which is defined 
as the time when the region within a half-mass radius 
achieves a two-body relaxation on average. Both the lin- 
ear TM relation and the LV relation appear even when 
the two-body relaxation is not attained. Hence, we can 
say that our two hypotheses are not affected through the 
two-body relaxation but are characteristic even in colli- 
sionless system in our estimation with t* el . 

From a viewpoint of these two hypotheses, we com- 
pared our results with two isotropic static models, i.e., 
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FIG. 13: The same as FigEDbut for King models with the 
unit ol G = th = Mtot = 1. The solid curve refers to the 
density for King model with c = 0.67. The other curves refer 
to the density for King models with c = 2.12 (dashed), c = 
2.72 (dotted). The plots represent the numerical data of SC 
(b t = 0,a = 0) with N = 10 4 . 
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FIG. 14: Deviation 5„ and St versus model parameters for 
polytropes and King models. Residuals with respect to poly- 
trope index n for polytropes and concentration parameter c 
for King models are plotted. Each curve refers to S v for poly- 
tropes (solid), 6 V for King (dashed), St for polytropes (dotted), 
St for King (dot-dashed). The gray zone refers to the range of 
Sv = 0.106 ± 0.053 derived from the numerical simulations of 
cold collapse in Fig^and Fig[3] Thick horizontal line denotes 
5^= 0.106. 



polytropes and King models. The polytrope with n = 5 
is special among those models, since it exactly satisfies 
the LV relation. Plummer's model, the polytrope with 
n = 5 with a flat boundary condition at the center also 
satisfies the constant TM relation quite well in the inter- 
mediate region 0.2M to t < M r < 0.8M tot . Moreover, its 
density fits well to the numerical data quite well in this 
intermediate region. On the other hand, both its TM re- 
lation and density profile in the outer region are slightly 



shifted from those in cold collapse simulations: the tem- 
perature falls off more steeply and the density behaves 
as r~ 5 in the outer region. As for King models, in con- 
trast, there is no parameter range in which both the TM 
relation and density profile fit well to the numerical data 
even in the intermediate region. 

We also examined the deviations from the two hy- 
potheses quantitatively for these two models. Only the 
polytropes with n ~ 5 are acceptable in that their devi- 
ations from both of the hypotheses are compatible with 
the numerical results of cold collapse. The LV relation 
imposes a strong constraint for these models. Especially 
for King models, the deviation from LV relation takes the 
maximum value around the parameter c ~ 2, in which 
the surface brightness fits well to i? 1 / 4 law. This sug- 
gests that considering only the asymptotic behavior of 
the bound state is not enough to describe the collision- 
less static models to express observational or numerical 
results. 

Combining the results for the above two models, we 
can speculate that Plummer's model follows both of the 
hypotheses we proposed (a) LV relation and (b) linear 
TM relation. However, this model is lack of the following 
properties of the bound states in the actual simulations 
of cold collapse: (i) deviation from a hydrostatic equi- 
librium state and (ii) anisotropy of velocity dispersion in 
the outer region. 

As for (i), the hydrostatic equilibrium condition is not 
attained especially when the violent degree of the ini- 
tial collapse is strong, since the bound state slowly ex- 
pands outward and is not balanced with a gravitational 
force after a cold collapse. Hence, to express these condi- 
tions we may have to consider not the static solution but 
the stationary solution admitting the coherent motion of 
matters. 

As for (ii), we derived that for a mild collapsing case, 
the LV relation fluctuates more largely in the outer region 
and the slope of the TM relation at this region becomes 
a little steeper at least in the collisionless stage. We sus- 
pect that these results are related to the anisotropy of 
this system. Hence for more precise description of the 
system possessing two hypotheses, it seems necessary to 
treat the model with the anisotropy of the velocity dis- 
persion. So far several collisionless static models with the 
anisotropic velocity dispersions or with the non-spherical 
symmetry have been proposed 0, E3, 123, E3, [S3, 
Incompleteness of the violent relaxation or the existence 
of conservative quantities during the phase mixing seem 
to be a candidate as the model to induce such anisotropy 
P, |23 • The hypotheses we proposed are expected to be 
used as a guideline to build up the models following the 
characters of the bound states after a cold collapse by 
combining with these anisotropic models. 

It also seems important how widely our two hypotheses 
are applicable in more general mixing processes including 
a merging process which yields a cuspy density profile in 
the central region |2f|. If we want to know the detail 
of such a tiny central region, such as the existence of 
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the cuspy density profile, we have to treat such a central 
region more carefully We speculate that in such a region, 
at least one of the two hypotheses we proposed will be 
broken. These points are now under investigation. 
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APPENDIX A: MASS FRACTION OF THE 
CENTRALLY FLAT REGION AGAINST THE 
OUTER ANALYTICAL SOLUTION © 

The analytical solution © has an inner unphysical re- 
gion at r < . Since the density goes to zero toward the 
center in this region, the cumulative mass in this region 
becomes very small. Here we will check how much the 
value is modified when we put the constant mass density 
in this region and connect it to the analytical solution 
at r = r*. From the continuity at r — r*, the density 
becomes 



P = 



Po (4/e) 



(r < r*) 



Po{ro/r) 4 e- r °/ r (r > r,) 



The cumulative mass at r — r* is 
The cumulative mass outside r* is 

4 



(Al) 



(A2) 



M out = 4tt / p 

/ro/4 V r 



( V) e -r /r r 2 dr = 47rpor 3 ^ _ g-4) 



Hence the mass fraction of the inner region becomes 

M in M in 4 



M in + M n 



3e 4 + 1 



(A3) 

s 

0.0243. (A4) 



APPENDIX B: DERIVATION OF THE LV 
RELATION FOR POLYTROPES 



Following JT^, we use the notation using the relative 
energy e and relative potential <j> defined as 



E := + 
(t> := + 



(Bl) 



respectively. For polytropes, using these notations, the 
distribution function is described as 



f(E) = 



Kl E n - 3 / 2 , (E > 0); 
0, (E < 0), 



(B2) 



For (IB2fl . density p and temperature T are simply ex- 
pressed as 



p = c n (p , 

, 3/cbT 
< v >= 



n + 1 



(B3) 
(B4) 



where c n := (2tt) 3 / 2 

From JSJl and {El}, we get 



W r = 



From and (|B4|l . we get 



^+db^ = <^TT) 



— 2tt / (j)p(r)r 2 dr. 



(B5) 



(B6) 



Especially for n = 5, from the condition $* = 0, the LV 
relation 2K r /\W r \ = 1 is obtained. 



APPENDIX C: DERIVATION OF TM RELATION 
FOR PLUMMER'S MODEL 

Following [l7j . we define the dimensionless variables 



ip := <f>/<f>0 s:=r/b, 



(CI) 



where (f>o is the value of the relative potential at the center 
and b := {AnGcn^ 1 } ~ 1 / 2 . Then, tp is followed by Lane- 
Emden equation, 



1 d 



,dib 



s 2 ds \ ds 
with the boundary condition 

lim ip(s) = ip 



-r, {i> > 0); 
0, (^ < 0), 



(C2) 



(C3) 



where ip* := $»/<^o- Especially for n = 5, the solution 
of (|C2|) with the flat boundary condition at the center 
becomes 



1 



and meets the boundary condition 

lim ij)(s) = V^* = 0. 

s — >oo 

The density is described as 

1 



(l + S 2 /3)"/2 



(C4) 



(C5) 



(C6) 
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In spherically symmetric case, the cumulative mass M r 
is expressed as M r = ^ whose dimensionless form 
becomes 



M r = -s 



ds ' 



(C7) 



and the total mass M tot is equal to Substituting (|C4I 



into dimensionless form of (|B4(1 and (|C7|I and eliminating 
s, we can easily derive the TM relation of Plummer's 
model as 



m 6 V \M tot 



2/3 



(C8) 
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